Breakdown of correspondence in chaotic systems: 
Ehrenfest versus localization times 
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Breakdown of quantum-classical correspondence is studied on an experimentally realizable ex- 
ample of one-dimensional periodically driven system. Two relevant time scales are identified in this 
system: the short Ehrenfest time tn ~ ln(7i _1 ) and the typically much longer localization time scale 
Tl ~ h~ 2 . It is shown that surprisingly weak modification of the Hamiltonian may eliminate the 
more dramatic symptoms of localization without effecting the more subtle but ubiquitous and rapid 
loss of correspondence at t%. 
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Quantum realizations of classically chaotic systems be- 
have differently than their fully classical counterparts. 
The breakdown of quantum-classical correspondence oc- 
curs on two distinct time scales: the short Ehrenfest time 
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gives the time scale at which the quantum minimal 
wavepacket spreads sufficiently over a macroscopic part 
of the phase-space to feel nonlinearities in the potential 
. In (|l|) Ap denotes the initial uncertainty in mo- 
mentum, A is the Lyapunov exponent, x — y/d x V/d^V 
is a typical scale on which nonlinearities in the poten- 
tial V are significant. Ehrenfest time tn marks the mo- 
ment when quantum corrections become important in the 
time evolution of the Wigner function W(x,p), governed 
by the Wigner transform of the von Neumann bracket, 
known as the Moyal bracket B : 



dt 



W = {H, W} MB 
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The Moyal bracket appearing above is defined as 



{H, W} MB 



The arrows indicate the direction of differentiation. 

In a chaotic system initially localized distribution will 
be stretched in the unstable directions. In a Hamilto- 
nian system, this stretching must be matched by shrink- 
ing in the stable directions as mandated by the Liouville 
theorem. Consequently, the typical width of the effec- 
tive support of classical probability distribution shrinks 
as exp(— Xt). That leads to an exponential increase of 



higher order derivatives of W in (||), so after t% quan- 
tum corrections become comparable to terms given by 
the Poisson bracket || . 

There is evidence that at t% also the averages over the 
phase space show symptoms of correspondence loss [||^] 
although not everyone agrees (7]|| . The first goal of our 
paper is then to investigate the dependence of the onset 
of the loss of correspondence in the averages on ft, and 
to show that it is consistent with (Q). 

In addition to the rapid correspondence loss on the 
logarithmic time scale t% there is plentiful evidence for 
the dynamical localization in a class of quantum chaotic 
systems |i|dTl|]. It sets in on a much longer localization 
time scale: 



TV = 



(3) 



where d is diffusion constant characterizing initial growth 
of variance in momentum of a quantum system. Tuto- 
rial example of localization is described in |12|. It is of 
obvious interest to consider systems in which both mech- 
anisms of correspondence loss are present, and to investi- 
gate their behavior with varying parameters. The system 
we investigate is a straightforward modification of ex- 
perimentally studied system based on an optical lattice. 
Therefore, it should be simple to verify our predictions 
in the laboratory. 

We study quantum-classical correspondence in a sys- 
tem described by the Hamiltonian 



p 

H(x. p) = k cos(x 



I smt) + 



(4) 



For a — this Hamiltonian reduces to a model recently 
investigated both theoretically |13-^] and experimen- 
tally pd| ] corresponding to motion of cold atoms in a 
quasi-resonant standing wave with periodically moving 
nodal pattern. Classical behavior of the system is of 
the mixed type with large regions of phase space domi- 
nated by a chaotic motion (e.g. Lyapunov exponent of 
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the chaotic component is A = 0.2 for chosen parameter 
values m — 1, k = 0.36, 1 — 3, a = 0, [H). 

To investigate quantum-classical correspondence one 
can start by looking at the Wigner function. Structures 
which appear in W are known to be different from (but 
not unrelated to) the classical probability distribution in 
the phase space. One can especially note existence of the 
interference fringes that saturate on small scales associ- 
ated with action s ~ (2irh) 2 /S, where S is the action of 
the system - e.g., the phase space volume over which W 
can spread |Q . 

One might nevertheless suppose that the quantum 
structures appearing on such small in the Wigner func- 
tion will have little effect on the averages. In particular, 
averages typically concern large scales, and, therefore, 
should not be directly effected by small-scale interfer- 
ence. This is, however, only partially true. 
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FIG. 1. Difference between classical, 



A 2 



and quantum, 



A qm , variances for three values of Ti and for the same initial 
condition i.e. coherent state centered at (0.0) in the phase 
space. In the insert, in dependence on log(l/ft) averaged over 
few initial conditions is shown. Error bars indicate standard 
deviation of the results for given value of Ti. Criterion used 
to estimate Ehrenfest time scale compares magnitude of Pois- 
son bracket and a quantum correction (the rest of the Moyal 
bracket in (g)). 

Figure [l] shows the difference between classical, A^, 
and quantum, Aq m , averages of a typical quantity (A 2 = 
(p 2 ) - (p) 2 ) for h = 0.16, 0.064 and 0.0256 obtained in 
the time evolution resulting from Hamiltonian (|J) with 
a = 0. Initial classical distribution is a two-dimensional 
Gaussian in phase space and corresponds to a Gaussian 
wavepacket used for quantum evolution. The loss of cor- 
respondence is apparent: early on, a noticeable difference 
between the two averages develops. The smaller the value 
of h , the later it appears. The insert shows that the 
discrepancies occur after a time which seems indeed log- 
arithmic in h. They persist without a significant change 



in their magnitude until localization (which occurs at a 
Ti, ~ h~ 2 , a time scale much more sensitive to changes of 
h than the logarithmic t%) begins to dominate. Vertical 
error bars in the insert arise from the fact that the local 
Lyapunov exponents (i.e. evaluated over a time scale of 
order t%, (0)), differ for various initial conditions at given 
h. This in turn affects value of t% (see Fig. ||). 




FIG. 2. Difference between classical and quantum vari- 
ances for h — 0.0256 and four initial conditions [coherent 
states centered at points (0.0), (vr/2, 0), (— tt/2, 0) and (n, 0) 
in the phase space] is plotted using black, red, green and blue 
lines respectively. The insert shows the discrepancy between 
classical and quantum averages as defined by (^) as a function 
of h. The dependence seems to be logarithmic. 

It is intriguing to enquire about the magnitude of the 
typical discrepancy between the quantum and classical 
averages in the time interval between t% and Tj,. To ad- 
dress this question we have evaluated 
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for tn < ti < t-2 < Tl as a function of Ti. The results 
(averaged over several initial conditions) show that 5(h) 
decreases with decreasing h significantly slower than lin- 
early (compare the insert in Fig ||). While our data do 
not allow for a definite conclusion, they point out towards 
a logarithmic ti dependence of 5(h). The results obtained 
for 5(h) depend slightly on the choice of t\ and t 2 - That, 
together with the size of the error bars, does not allow 
us to be more definitive. The logarithmic dependence is 
not in disagreement with results of p"7t (who claim that 
largest deviations are 0(h)) obtained for a model of two 
interacting spins. Our data indicate a similar behavior 
for a realistic, experimentally accessible system. Still, the 
latest result for the interacting spin system [jllj gives a 
square root rather than the logarithmic behavior. 

To the best of our knowledge, tn has not been yet de- 
tected experimentally. On the other hand, it is known 
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that random external noise 
E 



201 as well as decoher- 



22 1 suppress evidence of the breakdown at 



t% or at Tl, restoring quantum-classical correspondence. 
This happens when the coherence length [p2|j23| 



l c = Hy/X/2D 
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maintained by the chaotic system in the presence of 
decoherence (which can be parametrized under certain 
conditions J2^,|2j| by the addition of the diffusive term 
~ DdpW to the evolution equation, (Q)) becomes small 
compared to the scale of nonlinearities \ 



(7) 



The price for the restoration of quantum-classical corre- 
spondence by decoherence is irreversibility p,^2[. 

Let us now consider the localization. At the onset of 
localization (clearly visible in Fig. || for a = 0) the differ- 
ence between the classical and quantum variances begins 
to rapidly increase, in a manner approximately consis- 



tent with — Aq m = d(t - Tl). This is to be expected. 



In systems effected by localization wave functions ex- 
hibit characteristic exponential form, as exemplified in 
Fig. H(a). 




FIG. 4. Time average probability density in momentum 
representation: (a) exponentially localized for a — 0, (b) 
Gaussian for a — 0.01. 



In the presence of even a weak harmonic force a promi- 
nent growth of the difference between quantum and clas- 
sical averages nearly disappears. Also wave functions 
cease to be exponentially localized - compare Fig. ||(b). 
At the first blush, one may find this dramatic response 
to even a weak harmonic potential quite surprising. 
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FIG. 3. Difference between classical and quantum variance 
for four values of parameter a, as indicated in the figure, and 
for h = 0.16 . Addition of harmonic potential destroys dy- 
namical localization present for a = 0. All initial conditions 
are the same. 

Localization can be suppressed by noise and decoher- 
ence, as was recently shown by two groups using similar 
experimental setups |^3| . Figure |^ shows that a much 
simpler mechanism for suppression of localization exists: 
When the driven periodic potential is supplemented by 
even a weak harmonic potential, [controlled by the pa- 
rameter a in (Q)] localization disappears while chaotic 
character of motion persists, as illustrated by Poincare 
surfaces of section in Fig. ||. 




position 

FIG. 5. Stroboscopic Poincare surface of sections for the 
system with parameter a = (left), and for a — 0.01 (right). 

Our explanation of this effect is simple and general. 
Numerical studies show that the suppression of localiza- 
tion occurs when quarter of the period of the harmonic 
oscillator is approximately the localization time scale Tl- 
Localization can be maintained only when the "swap- 
ping" between x and p generated by the harmonic evo- 
lution happens on a time scale longer than Tj,. (After 
all, localization in p is incompatible with localization in 
x, since the wave function cannot have the exponential 
form in both complementary observables.) In fact, we 
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have seen disappearance of localization when the period 
of the harmonic oscillator is about WTl- So, the "swap- 
ping" does not need to be very frequent. This simple 
explanation of suppression of localization is model inde- 
pendent and may be the key result of our paper. 

Experimental verification of our two key results con- 
cerning (i) tfr and (ii) Tl should be possible. It will be 
more difficult to check for the intriguing prediction of the 
correspondence breakdown on short time scale tfr. One 
would ideally want to have otherwise identical two sys- 
tems, one of them quantum, the other classical. In our 
quantum Universe this is unfortunately impossible. One 
remedy would be then to compare real experimental re- 
sults with the classical computer experiment for the same 
values of all relevant parameters. 

There is, however, something disingenuous about com- 
paring a classical computer simulation with a quantum 
experiment in an attempt to find a breakdown of quan- 
tum classical correspondence. We shall therefore sug- 
gest an alternative strategy: Effective value of Ti can be 
adjusted in the experiments Jll],|23|]. Therefore, it may 
be possible to look for differences between some average 
quantities, say A2, for two otherwise identical systems 
which are quantum to a different extent that is adjusted 
by selecting different values of h. It is obvious from 
Fig. [l] that significant differences between two systems 
endowed with different effective value h would appear, 
and it should be possible to check for the validity of (Q). 

Destruction of localization by introduction of the har- 
monic trap into experiments which test quantum chaos 
in optical lattices should be much easier to verify. It 
may be desirable to superpose optical lattice of the sort 
used to investigate quantum chaos |IT| , p3| ] in a magneto- 
optical trap anyway (e.g. to confine the system such as 
BEC, where investigation of chaos may be of interest in 
its own right Hj). It should be then relatively easy to 
see if our prediction of the loss of localization is indeed 
borne out by the experiments. 

Both of our key conclusions, while investigated numer- 
ically in a specific system selected for its relevance to 
experiment should have a broad validity. In particular, 
destruction of localization due to rotation of phase space 
as a whole by the background (e.g. harmonic) potential 
on a time scale comparable to the localization time Tl 
is justified by a simple physical argument. This effect is 
dramatic and should be easily accessible to experiments. 
The earlier and more subtle phenomena that occur at tfr 
are robust to the addition of an overall potential, provid- 
ing that the chaotic character of the system is preserved. 
Our suggestion for investigation of the phenomena relies 
on the ability to adjust the relevant effective Ti in exper- 
iments. Moreover, we have found indication that such 
quantum-classical discrepancies in the averages diminish 
more slowly than 0(h), which - if confirmed in experi- 
ments we suggest - would be of major importance to our 
understanding of the classical limit of quantum theory. 



JZ and ZPK acknowledge support by KBN under grant 
2P03B00915. 



[f] G.P. Berman and G.M. Zaslavsky, Physica (Amsterdam) 

91A, 450 (1978). 
[2] G.M. Zaslavsky, Phys. Rep. 80, 157 (1981). 
[3] W.H. Zurek and J.P Paz, Phys. Rev. Lett. 72, 2508 

(1994); ibid. 75, 351 (1995). 
[4] M. Hillery, R.F. O'Connell, M.O. Scully, and E.P 

Wigner, Phys. Rep. 106, 121 (1984). 
[5] F. Haake, M. Kus, and R. Scharf, Z. Phys. B 65 381 

(1987). 

[6] S. Habib, K. Shizume, and W.H. Zurek, Phys. Rev. Lett. 

80, 4361 (1998). 
[7] W.A. Lin and L.E. Ballentine, Phys. Rev. Lett. 65 2927 

(1990). 

[8] R.F. Fox and T.C. Elston, Phys. Rev. E 49, 3683 (1994); 

ibid 50, 2553 (1994). 
[9] S. Fishman, D.R. Grempel, and R.E. Prange, Phys. Rev. 

Lett. 49, 509 (1982). 
[10] G. Casati, I. Guanieri, and D.L. Shepelyansky, IEEE J. 

Quantum Electron. 24, 1420 (1988). 
[11] F.L. Moore, J.C. Robinson, C. Bharucha, P.E. Williams, 

and M.G. Raizen, Phys. Rev. Lett. 73 2974 (1994); J.C. 

Robinson, C. Bharucha, F.L. Moore, R. Jahnke, G.A. 

Georgakis, Q. Niu, M.G. Raizen, and B. Sundaram, Phys. 

Rev. Lett. 74, 3963 (1995). 
[12] F. Haake, Quantum signatures of chaos, 2nd ed. 

(Springer, Heilderberg, 1992). 
[13] B.V. Chirikov and D.L. Shepelyansky, Zh. Tekh. Fiz. 52 

238 (1982). [Sov. Phys. Tech. Phys. 27, 156 (1982).] 
[14] R. Graham, M. Schlautmann, and P. Zoller, Phys. Rev. 

A, 45 R19 (1992). 
[15] P.J. Bardroff, I. Bialynicki-Birula, D.S. Krahmer, G. 

Kurizki, E. Mayr, P. Stifter, and W.P. Schleich, Phys. 

Rev. Lett. 74, 3959 (1995). 

[16] W.H. Zurek (in preparation). 

[17] J. Emerson and Ballentine, quant-ph/001102C . 
[18] J. Emerson and Ballentine, quant-ph/0103050. 
[19] E. Ott, T.M. Antonsen, Jr., and JD. Hanson, Phys. Rev. 

Lett. 53, 2187 (1984). 
[20] T. Dittrich and R. Graham, Ann. Phys. (N. Y.) 200, 363 

(1990). 

[21] P.A. Miller, S. Sarkar, and R. Zarum, Acta Physica 
Polonica 29, 3643 (1998). 

[22] W.H. Zurek, Physica Scripta T76, 186 (1998). 

[23] G.H. B all, K.M.D. Vant, H . Ammann, and N.L. Chris- 
tensen, |quant-ph/9902006| ; K. Vant, G. Ball, H. Am- 
mann, and N. Christensen, Phys. Rev. E, 59, 2846 
(1999). 

[24] H. Ammann, R. Gray, I. Shvarchuck, and N. Christensen, 

Phys. Rev. Lett. 80, 4111 (1998). 
[25] B.G. Klappauf, W.H. Oskay, DA. Steck, and M.G. 

Raizen, Phys. Rev. Lett. 81, 1203 (1998); Phys. Rev. 

Lett., 82, 241 (1998). 
[26] SA. Gardiner, D. Jaksch, R. Dum, J.I. Cirac, and P. 

Zoller, Phys. Rev. A 62, 023612 (2000). 



4 



